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FOREWORD 
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NOMENCLATURE 


A 

= 

cross-sectional area (in. 2 ) 

A t 

= 

flow area between seal knives and land, seal throat (in. 2 ) 

CL 

= 

clearance between seal knives and land (in.) 

DTC 

= 

distance to contact: axial clearance between knife and land, undefined for constant 
height straight-through seals (in.) 

f() 

= 

function of the variables ( ) 

f 

= 

Fanning friction factor 

gc 

= 

standard gravitational acceleration mass conversion factor (lb m ft/lb b ) 

H 

= 

seal height (in.) 

H 

= 

4A 

hydraulic diameter, H = (in.) 

Kc 

= 

contraction coefficient 

Ke 

= 

expansion coefficient 

Kf 

= 

wall friction loss coefficient 

KH 

_ 

knife height (in.) 

KN 

= 

number of knives 

KP 

= 

knife pitch (in.) 

KR 

= 

knife tip radius (in.) 

KT 

= 

knife tip thickness (in.) 

K vf 

= 

Venturi-friction coefficient 

Kp 

= 

knife taper angle (degree) 

K 0 

= 

knife slant angle (degree) 

L 

= 

length of seal (in.) 

LTSD 

= 

leakage flow direction from large-to-small seal diameter 

M 

= 

mach number 

P 

= 

wetted perimeter of duct (in.) 

Ps 

= 

local static pressure (psia) 

Pd 

= 

seal plenum downstream pressure (psia) 

Pr 

= 

seal pressure ratio, Pu/Pd 

Pt 

= 

local total pressure (psia) 
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r 

R 


Re 

SH 
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Tu 

V 

w 

XMUL 

oc 

T 

5 

e 

F 

K 

P 

♦ 


= seal plenum upstream pressure (psia) 
= Pd/Pu 


= gas constant 


r lb f ft ) 


= streamwise Reynolds number, 


pVH 

F 


= step height (in.) 

= leakage flow direction form the small-to-large seal diameter 
= local total temperature (°F) 

= seal upstream plenum temperature (°R) 

= leakage gas velocity (ft/sec) 

= seal airflow rate (lb m /sec) 

= area correction factor for clearance above a knife downstream of a step 

= jet expansion angle (degree) 

= ratio of specific heats 

= jet expansion height (in.) 

= land surface roughness (pin.) 

^ J . . . ( lb m 

= fluid dynamic viscosity — — 

V ft sec 

= conventional transcendental number, ratio of circular circumference to diameter 


= density 



PuA, 


airflow parameter 


lb. 


lb 


O R 1/2 ^ 

f sec j 
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1.0 INTRODUCTION 


NASA's advanced engine programs are aimed at progressively higher efficiencies, greater 
reliability, and longer life. Recent studies have indicated that significant engine performance 
advantages can be achieved by employing advanced seals [1]*, and dramatic life extensions can 
also be achieved. Advanced seals are not only required to control leakage, but are necessary to 
control lubricant and coolant flow, prevent entrance of contamination, inhibit the mixture of 
incompatible fluids, and assist in the control of rotor response. 

Recognizing the importance and need of advanced seals, NASA, in 1990, embarked on a five-year 
program (Contract NAS3-25644) to provide the U.S. aerospace industry with computer codes 
that would facilitate configuration selection and the design and application of advanced seals. 

The program included four principal activities: 

1 . Development of a scientific code called SCISEAL, which is. a Computational Fluid 
Dynamics (CFD) code capable of producing full three-dimensional flow field 
information for a variety of cylindrical configurations. The code is used to enhance 
understanding of flow phenomena and mechanisms, to predict performance of complex 
situations, and to furnish accuracy standards for the industrial codes. The SCISEAL 
code also has the unique capability to produce stiffness and damping coefficients that 
are necessary for rotordynamic computations. 

2. Generation of industrial codes for expeditious analysis, design, and optimization of 
turbomachinery seals. The industrial codes consist of a series of separate stand-alone 
codes that were integrated by a Knowledge-Based System (KBS). 

3. Production of a KBS that couples the industrial codes with a user friendly Graphical 
User Interface (GUI) that can in the future be integrated with an expert system to 
assist in seal selection and data interpretation and provide design guidance. 

4. Technology transfer via four multiday workshops at NASA facilities where the results 
of the program were presented and information exchanged among suppliers and users 
of advanced seals. A Peer Panel also met at the workshops to provide guidance and 
suggestions to the program. 

This final report has been divided into separate volumes, as follows: 

Volume 1: Executive Summary and Description of Knowledge-Based System 

Volume 2: Description of Gas Seal Codes GCYLT and GFACE 

Volume 3: Description of Spiral-Groove Codes SPIRALG and SPIRALI 

Volume 4: Description of Incompressible Seal Codes ICYL and IF ACE 

Volume 5: Description of Seal Dynamics Code DYSEAL and Labyrinth Seal Code KTK 

Volume 6: Description of Scientific CFD Code SCISEAL. 


Numbers in brackets designate references presented in Section 7.0. 
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This volume describes two codes: DYSEAL and KTK. The code DYSEAE determines the 
response of face and cylindrical seals to shaft excursions. The code KTK is a labyrinth code that 
was obtained from the Air Force and produces leakage information for straight and stepped 
labyrinth seals and includes routines for geometric optimization. References 2 and 3 provide the 
details of the code implementation. Reference 3 is an edited version of Volume IV of report 
AFWAL-TR-85-2103 prepared by Allison Engine Company for the U.S. Air Force. 

1.1 Code DYSEAL (Dynamic Response of Seals) 

Dynamic response of seal rings to rotor motions is an important consideration in seal design. For 
contact seals, dynamic motions can impose significant increases in interfacial forces, resulting in 
high wear and reduction in useful life. For fluid film seals, the rotor excursions are generally 
greater than the film thickness, and if the seal ring does not track, contact and failure may occur. 
The computer code DYSEAL can determine the tracking capability of fluid film seals and can be 
used for parametric geometric variations to find acceptable configurations. 

The type of seals that can be analyzed are depicted in Figures 1 through 3 . Figure 1 shows a 
stationary seal ring and a rotating mating ring. The secondary seal is a piston ring with radial 
pressure loading on the OD. The shaft or rotor can be given five degrees of freedom, consisting 
of three translations (x, y, and z) and two rotations about the x and y axes, respectively. The seal 
ring response is also in five degrees of freedom. The interface is represented by cross coupled 
stiffness and damping coefficients that are obtained from other codes. The effects of Coulomb 
friction of the secondary seals on seal ring response are included. Figure 2 shows an inverted 
configuration with the initial radial pressure on the piston ring on its ID. This inside configuration 
results in less pressure loading on the ring because the ID area is less than the OD area. The 
reduced loading also reduces the secondary seal ring friction that may retard tracking. In addition 
to piston ring secondary seals, an O-ring secondary can also be applied. 

Figure 3 shows a floating ring seal that can also be analyzed by the code. This configuration 
permits two degrees of freedom for both the shaft and ring, and is intended to determine seal ring 
response to an orbiting shaft. The secondary seal occurs between the ring and the wall and x-y 
Coulomb friction at that location is accounted for. 

The method of computation is a forward integration in time that provides absolute motions in all 
degrees of freedom. The reason that this approach was chosen was because of complications 
caused by Coulomb friction. At every time step, friction has to be evaluated to determine if 
motions continue or are halted. 

1.2 Code KTK 

The computer code KTK calculates the leakage and pressure distribution through labyrinth seals based 
on a detailed knife-to-knife (KTK) analysis. Input data are required to describe in detail the seal 
geometry and the environmental conditions affecting the leakage. Output is provided in the form of 
leakage flow and flow resistance characteristics, i.e., flow factor versus pressure ratio. In addition, an 
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optimization feature is included which permits the user to identify global geometric constraints and 
allows the code to identify an optimum seal configuration based on minimum leakage. 

This volume describes the fundamental theory of the design model, features of the design model 
computer program KTK, and resulting output data. Appendix A includes sample output data. 



831834 -2 


Figure 1. Fluid Film Face Seal Parameters 
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2.0 THEORETICAL DEVELOPMENT OF CODE DYSEAL 


The code determines the response of the seal ring in five degrees of freedom to shaft vibrations in 
as many as five degrees of freedom. These degrees of freedom are: 

1 . x & = seal ring displacement in x direction. 

2. y s = seal ring displacement in y direction. 

3 . z s = seal ring displacement in z direction. 

4. p s = seal ring rotation about x-x axis. 

5. a & = seal ring rotation about y-y axis. 

Note that throughout this report, seal motions are subscripted with an s and shaft motions are 
nonsubscripted. Unit vectors are j ? and k in the x, y, and z directions, respectively. Coulomb 
friction is accounted for in both the secondary seal and the interface. 

2.1 Equations of Motion 

Considering small motions, the following equations apply: 


SF^ = mx s 

(2.1) 

XF y = my. 

(2.2) 

XL = mz, 

(2.3) 


(2.4) 

EM y = I,a, 

(2.5) 


where m = mass of seal ring and I — transverse moment of inertia of seal ring. 

2.2 Development of Newmarks’ Method 


The solution to the equations of motion are obtained by the use of Newmarks 1 method or the 
average acceleration method [4], The velocity u i+1 at a time station, i+1 is approximated as 




f u,+0 M ' 


At 


Similarly, 


u lt ,=u,+ 


A.+U,,, 


At 


(2.6) 


(2.7) 
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If the value of u from Equation (2.7) is substituted into Equation (2.6), the following equation 

results: 


U i+1 =U i +U i At + 




At" 


( 2 . 8 ) 


but 


U i+1 =M- l {F 1+ ,-KU ltl -DU it ,} 


(2.9) 


where M = mass matrix; F = unbalanced force vector; K = stiffness matrix; and D = damping 
matrix. 


Substituting Equation (2.9) into (2.8) produces: 


(■ — + KM" 1 ) = ErU, +— 0, + U. + M"'F i+l - DMT'U 


2 ■ ■ At 


i+1 


Ui+ W ) At 

Now from Equation (2.7), 

tt = 2U i±L _2U L _ T -T 

i+1 At At 1 

Substituting Equation (2. 1 1) into (2.10) and multiplying by M produces: 


4M 

At 2 


2D\ 


( 4M 2D V ( 4M 


+ K + ^- J U i+i - F i + i+l — + — l U ! + 


At At ) 


At 


+ D U ; +MU; 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


Thus, an expression has been derived that relates the displacement at the new time step to 
displacements, velocities, and accelerations at the prior time step. 

Once u ;+1 is obtained, u and u i4 ., are obtained from Equations (2.7) and (2.6), respectively. 




( 2 - 13 ) 


From Equation (2.6), 


u i+1 =^[u i+ , -U,]-U i 


(2.14) 
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Substituting Equation (2.13) into (2.14) gives: 


4 / s 4 • 

U‘.. = -(u i+1 - U,)-— U i -U i 


At' 


At 


(2.15) 


Thus, displacements, velocities, and accelerations are determined from the results of previous time 
steps. Initially, these quantities equal zero ultimately. 


2.3 Solution Process 


Figure 4 is a flow diagram of the program logic. The program computes the mass and inertia 
properties of the seal ring and the location of the center of gravity. After computing all constants 
and matrix elements that are independent of time, the program enters the time step loop. Shaft 
motions are incremented first. Using updated shaft motions, the secondary seal friction is 
determined. This includes friction magnitudes and direction in the x , y s „ z , (3 s , and a § directions 
as well as the friction components that go into the stiffness and damping matrices and force 
vector. 


The force vector, F, is next updated because, as indicated in using Newmarks method, the most 
recent force vector, F. is required. Then, Newmarks method is applied and the new seal 
displacements, velocities, and accelerations are determined. Subsequent to the calculations, 
adjustments are made to these variables because of friction resistance. The following paragraphs 
describe the development of the theory for the individual steps in the solution process, as outlined 
in Figure 4. 

2.4 Initialization 

Displacements, velocities, and acceleration are initialized prior to entering the time step loop. 
Initial displacements correspond to the shaft displacements at the first time step, so that the seal 
ring and shaft are in correspondence. Initial values of velocities and acceleration are nulled. 

2.5 Mass Matrix 


The code develops the mass and inertia properties from a series of connected ring elements. Up 
to 20 elements can be inputted with individual OD, ID, length, and density. From this input, the 
code determines the location of the center of gravity (CG), the mass, and the polar and transverse 
moments of inertia of the seal ring. Computed values are included in program output. 
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2.6 Computation of Constants 

The closing area is the unbalanced hydraulic closing area that varies for the type of seal being 
analyzed. For a piston ring seal (Figure 1), the closing area is: 

A C l = 7c ( R o»- R L) ( 2 .16) 

If an inside ring is employed (see Figure 2), then 

A CL = n(K-Rl) (2.17) 

For an O-ring secondary seal, there is no distinction between the inside and outside radii and the 
closing area is given by: 

A cl =*( r L- r L) (2' 18 > 


The same expression applies for a ring seal except that R sc is taken as the inside radius of the seal. 
The interface area is the mating area and is given (for all seals) by: 

A 1f =*( r L- r ? s ) (2' 19 ) 

Another area of interest is the difference between the interface and closing areas. The absolute 
difference between the interface and closing areas is: 


^CLL = I^CL — A IF I 

For a face seal, the hydraulic closing force is: 

F -PA + p a 

1 HCL z H^CL ^ r L^CLL 

where P H = high pressure and P L = low pressure. 
For ring seals. 


( 2 . 20 ) 


( 2 . 21 ) 


FHCL = R HA a .-P L A CLL (2-22) 

The code computes the spring preload by summing the preload from the individual springs. For a 
piston ring, there will be initial preloads from the installation spring stiffness and from pressure 
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on the ring circumference and on the ring face. The secondary seal preload per unit length for a 
ring pressurized on its OD is: 


p re( =p H ^+(pJ, 


(2.23) 


where p ^ = preload per unit of circumference; R q = outside radius of the piston ring; R = inside 
radius of the piston ring; w - the width of the piston ring; and (p \ - installed preload per unit 
of circumference. 


The total preload is: 

P r = 27tR i P re< (2.24) 

For an inside ring, R L and R q are reversed in computing p . 

For an O-ring seal, the preload per unit of circumference is an input quantity and the preload is 
given by Equation (2.24). 

The face seal axial forces are a function of preload and the coefficient of friction, such that 

Ff Z -P r xv (2.25) 

where F fe = secondary seal friction in the axial direction and v = coefficient of friction. 

The initial interface preload includes components from closing pressure and spring load and is 
equal to: 

Fif = F hcl + F sp (2.26) 

where F ff = initial interface load and F sp = initial spring closing load. 

The code computes the initial axial position of the seal ring accounting for secondary seal ring 
friction. For face seals, the equilibrium fluid film interface force is an input quantity. The 
procedure is to balance the closing loads by the fluid film load using the axial stiffness of the fluid 
film to determine position iterations. Figure 5 shows the algorithms used. 

For a ring seal, the fluid film stiffness is replaced by the structural stiffness of the seal ring. Often, 
a soft material such as carbon is used for the seal ring and its initial compression and face load are 
of interest. The ring seal stiffness is approximated by: 

K ZZ =A 1 F E/L 

where A JF = interface area; E = elastic modulus of the seal ring; and L = seal ring length. 
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2.7 Stiffness and Damping Outside of the Time Step Loop 

There are stiffness and damping quantities that are invariant and can be matrixed outside of the 
time step loop. 

2.7.1 Fluid Film Stiffness and Damping 

The fluid film interfaces are represented by cross coupled stiffness and damping coefficients 
that are obtained from other codes. For face seals, the fluid film has three degrees of 
freedom (z, (3, a) and the stiffness and damping quantities occupy the lower right portion of 
the 5x5 stiffness and damping matrix. The ring seal fluid film has two degrees of freedom 
(x and y) and the stiffness and damping values occupy the upper left portion of the matrices. 
Tables 1 and 2 show the locations of the stiffness and damping quantities. 

2.7.2 Spring Stiffnesses 

The total spring force is (see Figure 6): 

Nsp 

F !p =k SP 5X (2.27) 

1=1 

where F sp = total spring force; k sp = spring stiffness; §^ p = displacement of ith spring; and 
N = number of springs. (Note: the spring preload does not enter into the equations of motion 
since it is equilibrated by initial conditions.) The displacement of the ith spring is: 

8^=(8,+?.x^)-k (2.28) 


where 


5 = displacement of eg = \ 


<() s = rotation of seal ring about axis through eg = 


fPl 

a 

o 


= vector from eg to ith spring = < 


R 


sp 


R 


sp 


cosG 1 
sin0 l > 


z 


sp 
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Table 1. Stiffness Coefficients 


F \ A X Y Z P 


F, 

K K 
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◄ — Ring Seals 
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IYy X ^yy 


F, 


K zP 

k 2 « 

M, 

Face Seals — ► K pz 

Kpu 


M y 

K« 


, K »« 
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Table 2. Damping Coefficients 


\ A 

X 

Y Z 

P 

a 

F, 

Fy 

D„ 

Dy* 

D yy 
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D yy 


Fy 


D, 

D z p 

D„ 

M, 

Face Seals — D pz 

D pp 

■Dpa 

M y 


D„ 

D a p 

D aa 


95TR34-V5 


NAS A/CR— 2004-2 1 3 1 99/VOL5 


11 



The moments of the spring forces about the center of gravity are: 


M. =y r ' x F ' = 

sp sp sp 


M 

M 


yy 


o 


The axial stiffness of the spring is: 


N • k 

sp sp 


Rotational stiffness can be obtained explicitly. 


Aj 


*3 = 


where i = moment axis and j = displacement. 

For a single spring, the rotational spring constant is: 


(2.29) 


(2.30) 


(2.31) 


K 


spt 


k R 2 

sp sp 

2 


(2.32) 


The program numerically computes the spring stiffnesses and then adds them to the stiffness 
matrix for Newmarks computations. 


Stiffness and damping are also computed for the O-ring secondary seals, and is presented in 
Section 2. 14 along with the discussion of O-ring friction, which is a parameter whose direction 
varies with time. 


2.8 Shaft Increments 


Shaft motions are incremented inside the time step loop according to the following equations: 


x — x 0 cosco x t (2-33) 

y = y 0 sin co y t (2.34) 

z = z o sin co z t (2.35) 

P = |3 o cosC0pt (2.36) 

a = a o sinco a t (2.37) 


where x = shaft displacement in x direction; y = shaft displacement in y direction; z = shaft 
displacement in z direction; P = shaft rotation about x axis; and a = shaft rotation about y axis. 
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To simulate circular orbits, x and y are 90° out of phase. The amplitudes x q , y etc., are input 
quantities and can be arbitrary to simulate elliptical shaft orbits. Also, the frequencies of 
vibration, co^, etc., are input quantities and can be varied arbitrarily at the discretion of the 
user. Velocities and accelerations are computed by taking derivatives in the usual manner. 

2.9 Updating [K] and [D] 

For a ring seal, the fluid film stiffness and damping are constant quantities but their 
components in x and y vary with the position of the shaft, and thus they must be updated 
inside of the time step loop. Basically, the input values of K , K , etc., are values that are 
parallel and normal to the eccentricity vector. Referring to Figure 7, the position of the 
eccentricity vector varies as the shaft orbits. As shown in Figure 7, the eccentricity is along 
the x' axis. Then, for the primed axes, 

F' = -K'b'-D'b' (2.38) 


where 



which are input quantities 




The forces along x', y' must be transposed along x and y 
F= AF' 


( 2 . 39 ) 
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where 


cos0 -sin0 
A = 

sin 0 cos 0 

L_ J 

Substituting Equation (2.38) into (2.39), we obtain 

F = A[-k'8'-D'8'] 

but 6' = A T 8 and 8 = A T 8. Therefore, 

F - -AK'A t S - AD'A T 8 
but, F also equals 

F = -KS-DS 

where 


k 

K 

XX 

xy 

K 

K„„ 

yx 

yy 

D„ 

D* 

Dyj 

D yy 


(2.40) 


(2.41) 


(2.42) 


Therefore, comparing Equations (2.41) and (2.42) 

K = AK'A t and D = AD'A T (2.43) 

The code determines the position of the eccentricity vector by calculating the position of the 
minimum film thickness. The stiffness and damping transformations are appropriately added to 
the stiffness and damping matrices for NEWMARK computations. 
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2.10 Viscous Shear Forces and Moments 

For face seals, viscous shear forces are produced at the interface. These forces are 


f^ zH av 

h 



where 


(2.44) 

(2.45) 


(2.46) 


A = interface area 
V = seal ring velocity in x-y plane 
h = film thickness 
p = absolute viscosity 

The coefficients pA/h are included in the damping matrix. 

2.11 Applied Forces 

The computation of applied forces and moments are necessary for subsequent friction 
computations. The applied force vector includes all forces and moments excluding equilibrium 
and friction forces and moments. 

For ring seals, 

F, = -K(8.-8)-d(8.-8) (2.47) 


where: 


F a = applied force vector = 


ay 


K = film stiffness matrix 


K, 

K. 


xy 


yy 
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D = film damping matrix - 


D xx D xy 

D vx D, 


5 S = seal ring displacement = \ “ \ 

IM 

rsj 

5 = shaft displacement = < x > 

l S yJ 

Similarly, ^ and § are seal ring and shaft velocity vectors, respectively. 
For face seals, the matrix formulation is: 


F = -K8 -D8. 


where: 


F a = applied force vector = < 


ay 


Fa 

M 

M 


ay J 


8 r = relative displacement vector = 


X s 

y s 

z s ~z 

p s -p 

a -a 


5 r = relative velocity vector = < 


X s 

y s 

p s -p 

d-d 


(2.48) 


The forces from x and y displacements occur between the secondary seal and housing and are not 
relative with respect to the shaft. 


NAS A/CR— 2004-2 1 3 1 99/VOL5 


16 



Subscript, s, refers to the seal ring. Displacements without subscripts refer to the shaft motions. 
The elements of the stiffness and damping matrices are: 

K u = — K“ — O-ring stiffness in x direction 
K 22 = _k^ = O-ring stiffness in y direction 

K33 = -K z p + K zz = spring stiffness + fluid film stiffness in z direction 
K 34 = K^p = cross coupled film stiffness 
K 35 = K za = cross coupled film stiffness 
K 43 = Kp z = cross coupled film stiffness 

K 44 - _Kp p + Kp p = spring rotational stiffness + film stiffness about x axis 
K 45 = K = cross coupled film stiffness 
K 53 - K = cross coupled film stiffness 
K_ 4 = K = cross coupled film stiffness 

K 55 = -K^ p + K aa = spring rotational stiffness about y axis + film rotational stiffness. 

The damping matrix includes the viscous shear damping: 

Du = d sh -D“ = shear damping coefficient + O-ring damping in x direction 

D15 = -D SH z C g = shear damping coefficient x axial distance to eg 

D 22 = D sh -D“ = shear damping coefficient + O-ring damping in y direction 

D 25 = D SH z cg = shear damping coefficient x axial distance to eg 

D 33 = = film damping coefficient 

D 34 = D z p = cross coupled film damping coefficient 

D 35 = D = cross coupled film damping coefficient 

D z42 = D SH z cg = shear damping coefficient x axial distance to eg 

D 43 = Dp z = cross coupled film damping coefficient 

D44 = Dp p +D SH zJ g = film damping coefficient + shear damping coefficient x 
the square of the eg distance 

D 45 = D„ + D SH z^ g = film damping coefficient + shear damping coefficient x 
the square of the eg distance 

D m = -D SH z cg = shear damping coefficient x the distance to the eg 
d 3 = d « = cross coupled film damping coefficient 
D .4 = = cross coupled film damping coefficient 

D 55 = D aa +D SH z c 2 g = cross coupled film damping + shear damping coefficient x 
the square of the eg distance. 
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2.12 Piston Ring Secondary Seal Friction Forces and Moments 

2.12.1 Friction Forces and Moments from the Radial Surface of the Piston Ring 

The piston ring moves with the shaft in x and y, and can also hold back the shaft from moving 
(see Figure 8). Surface 2 of the piston ring is the radial face, and Surface 1 is the interior 
cylindrical surface. The velocity of Surface 2 is: 


v S c 2 = V + y s j + eo s xr 2 

“ s =Pj + aj 

^ A. As 

r 2 = Zjk + r f cos 0i + r f sin 0j 

■ A A. A /V 

G) s x r 2 = -P s z t j + p s r f sin 0k + d s Zj i - a s r f cos 0k 
It is assumed that there is zero k velocity of the piston ring. Therefore: 

v „2 = (*. + a s z, )i + (y, - Mi) J = Vi^ + VL 2y 

The direction of the friction force is opposite to that of the velocity 

- V 

Fjcj = friction force = -Po A P v 17f 1 7 

| ’ sc2 1 

where: 


P 0 = applied pressure on piston ring 

A p = unbalanced contact area of piston ring 

v = coefficient of friction 


and 


F sc2x = -Po A p V(S ‘ gnVSCZX) 

F S c2 y = -Po A p V(S ‘ gnVSCZy) 

The moment about the CG from the face friction force is: 

M 2 = f 2 X F sc2 = (z,k + r f cosOi + r f sin 0j) x (F„ 2 ,i + F„ 2y j) 
= z 1 F sc2x j -z 1 F sc2y i + k components that are neglected. 

where F „ and F , are defined above. 

sc'2x sc2y 


(2.49) 

(2.50) 

(2.51) 

(2.52) 

(2.53) 


(2.54) 

(2.55) 


(2.56) 
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The friction forces and moments are subsequently added to the force vector in the Newmark 
formulations. 


2.13 Friction Forces From the ID of the Piston Ring 

At the ID piston ring interface, there is only velocity in the z direction, which equals the relative 
velocity of the seal ring in the z direction. The major contribution to the normal force at the ID is 
the pressure that p o applies to the OD. The following equation results for the load per unit length 
at the ID of the piston ring. 


p = Po R sco W + p/ 
R , 


(2.57) 


where: 


p = preload per unit length of ID of piston ring 

p o = pressure on OD of piston ring 

R = outside radius of piston ring 

R d = inside radius of piston ring 

w = width of contact surface at ID 

p^ = initial or installed preload per unit length 

The direction of the friction force is opposite the direction of the axial velocity, y . If y =0, the 
direction of the friction force is opposite the direction of the applied force in the z direction, p . 
Therefore, 

F f = -2 jcv R scj P rf (sign V z or sign F^ ) (2.58) 


2.14 O-Ring Secondary Seal Stiffness and Friction Forces and Moments 

An O-ring secondary seal contributes stiffness and damping and friction forces and moments. 
Explicit analysis was conducted to determine these contributions. The O-ring is divided into 72 
segments of 5 degrees each. The displacement of the fth segment is: 

5* =u + (j>xr s f e (2.59) 

where: 


K 


N 


(V 

fR cosG^ 

1 sc 

<K 

>, u = < 

Un 

ii 

-©- 

«s 

r ,' = I sin 0' • 

K 


l°J 


0 

( Z sc -U (5) 
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The normal vector at 0^ is n* (Figure 9), where 


(2.60) 

(2.61) 


rr = cos0 £ i + sin0*j 

5*-n* =8j 

where = normal displacement at 0* . Similarly, the velocity of the seal ring at the ^th segment 
is: 


6 £ = u + <j> x r/ c 
where (2.62) 


V 


P.' 

U 2 

S (t> = < 

<*s 

0 


0 


and 


(2.63) 


(2.64) 


= O-ring stiffness per unit length 
D e ^ = O-ring damping per unit length 
Pr = O-ring preload per unit length 
F* = normal O-ring force at 0^ 


8' n =5'-n* 

The normal force at the fth segment on the seal ring is: 

K = -MXde' -D^sXde' -Pg-R^de' 

where: 


and 


K =FnCose' 

(2.65) 

F;=F'sin0' 

(2.66) 

f x = £f' 

C=1 

(2.67) 

72 

F y =lF y ' 

(2.68) 


f=\ 
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where: 


F* = x force at 0^ 

x 

Fy = y force at 0^ 

F x = total O-ring x force 
F y = total O-ring y force 

The O-ring moment is: 

M‘=f'xF n ' 


(2.69) 


and 


M = < 


n 


M 

M y 

Ml 


where: 


M* = moment due to the normal force at 0* 


and 


M x =fX, M y = fx 

1 = 1 

The stiffness of the O-ring is: 



(2.70) 


(2.71) 


where: 


K„ = ith stiffness in the i direction due to a j displacement 
5. = displacement in j direction. 

These stiffnesses are determined from a zero displacement position and are constant values. They 
are added to the stiffness matrix that is used in the NEWMARK computations. Similarly, 


D 


ij 



(2.72) 
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where: 


D.. = damping in i direction due to a j velocity. 

The O-ring friction imposes additional forces and moments on the seal ring. The friction forces 
are always along the z axis and direction is always opposite the velocity vector. The relative 
velocity at sector £ is: 

V'=(a'4xr')k (2.73) 

Then the friction force at the ^th segment is: 

F'=vF'(-signV') (2.74) 

where v = the coefficient of friction. The total friction force is : 


72 

F f =X F ' 

e = i 

The moment due to the friction force is: 


M f = i(r'xF/) 

e=i 



O' 


[ M d 

N 

luT 

0 

M f =< 





f 

o 


These forces and moments are added to the 


(2.75) 

(2.76) 

(2.77) 

vector used in the NEWMARK computations. 


2.15 Computation of the Force Vector 

The force vector contains all terms that are not directly multiplied by the seal ring displacements, 
velocities, or accelerations. 

mx s + kx s + Dx s = F (2.78) 
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where: 


m = mass matrix 
k = stiffness matrix 
D = damping matrix 
F = force vector 

x , x , x s = seal ring displacement, velocity, and acceleration vectors 

The terms of the force vector, F, contain stiffness and damping multiplied by shaft displacements 
and velocities plus friction restraint force: 

F = kx + Dx + F f 

where 

F = friction vector 
K = fluid film stiffness matrix 

D = fluid film damping matrix plus viscous shear damping terms 
F = friction restraint forces in all degrees of freedom 

2.16 Friction Restraint 

After the displacements, velocities, and accelerations are updated, it is necessary to determine whether 
friction should have halted the motion. When detennining friction restraint, it is important to realize 
that total forces and velocities are applied, and considering components alone can be misleading. For 
example, a body moving in a plane will not be restrained in a component direction as long as the total 
applied force exceeds the total friction force even though a component friction force can exceed a 
component applied force. For purposes of illustration, a piston ring secondary seal will be discussed. 
Referring back. Figure 8 showed the piston ring model. 

Wall friction on the piston ring will restrain lateral (x and y) and angular (a and p) motions of the 
seal ring. The ID friction of the piston ring will restrain axial motions of the seal ring. Consider a 
velocity versus time plot, as shown on Figure 10. 

There are three regions of accountability: 

1 . When accelerating, F f and F a are opposite |F f | < |F a | 

• F = friction force 

• F a = applied force. 

2. When decelerating, F f and F a are of the same sign; a finite velocity implies |F t | < |F a | - 

3. If the velocity changes sign between successive time steps, then somewhere between 
motion has stopped and cannot restart until |F a | > |F f |- Thus, there is a discontinuity in 
the velocity curve. If we followed the normal procedure without taking into account 
the finite stopping time, the velocity would be repositioned to point B in Figure 10 
instead of point A. 
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At the piston ring wall, the velocities of the ring in the x and y directions are: 


V x = x + az! 

( 2.79) 

V y =y-pz, 

(2.80) 

and the total velocity is: 


V T = A /V' +V y 

(2.81) 


where: 


V x = x component of velocity of piston ring 
V y = y component of velocity of piston ring 
x = x component of velocity of seal ring at CG 
y = y component of velocity of seal ring at CG 
a = rotational velocity about y axis 
p = rotational velocity about x axis 
z = axial distance from CG to piston ring wall 

The forces that would move the piston ring along the wall (x-y plane) come from both the lateral 

applied forces at the CG and the total applied moments about the CG. 

F«,=F„+M )y /z 1 (2-82) 

F ayt = F >y — / Zj (2.83) 


where: 


= x component of applied force at CG 
F ay = y component of applied force at CG 
= applied moment about y axis 
= applied moment about x axis 
F = total x component of applied force at ring wall 
F ayt = total y component of applied force at ring wall 

The total applied force at the ring wall is F , defined as: 

F,=V F » +I V ( 2 - 84 > 

Other parameters include: 

V„ =V* 2 +y 2 (2.85) 

where V = total transiatory velocity at CG. 
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( 2 . 86 ) 


F axy - V F ax + F ay 

where F = total translatory applied force at CG. The friction force at the wall is defined as: 

Ff=V F fc + F £ ( 2 ' 87 > 

where: 

F = total friction force 

F. = friction force in x direction 

fx 

F = friction force in y direction. 

xy J 

With these terms and definitions, a flow chart of the friction wall restraining algorithm is indicated 
on Figure 11. Note that even though friction may not restrain piston ring motion, a check has to 
be made on x and y motions at the CG. This is because the piston ring can move due to angular 
rotations about the CG without x, y translations of the CG. 

A similar routine has been established for restraint in the z direction. Since this is a single-degree-of- 
freedom motion, it is a much simpler algorithm than for the coupled x, y and angular modes. 

2.17 Minimum Film Thickness 

At each time step, the minimum film thickness is computed to determine if a negative film occurs. 
If so, the computations are halted and the seal is considered failed. 

Face Seal. Because of angular rotations, the minimum film thickness is computed at the outside 
radius of the fluid-film interface. The film thickness varies around the circumference of the seal 
and is thus a function of theta. In computing the film thickness, the circumference is subdivided 
into 72 increments, and the film thickness determined at each incremental intersection. 

H p =H 0 +AZ + <j>xT p -k (2.88) 

where: 

H q = equilibrium film thickness 

A Z = difference between seal ring and shaft displacement 
(j) = rotation vector 

r - position vector from center of gravity to point p 
For a face seal with a stationary (nonrotating) seal ring, 

H = H 0 +(Z s -Z) + [((p s -p)i +(cx s -a)j)x(-z cg k + R 0 cos0 p i + R 0 sin9 p j) k 
H = H 0 +(Z S -Z)+[(p s -P)R 0 sin 0 p - (a, ~a)R 0 cos0 p ] (2.89) 
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Ring Seal Clearance. Ring seal parameters are shown on Figure 12. Ring seals are limited to 
two degrees of freedom, x and y. Since both the shaft and ring can move, the film thickness is a 
function of relative displacements between them. The equations are as follows: 

<l>, = <h-<l> 

H = C + ^ r n 

(2.90) 

where: 


= displacement vector of seal ring 
£ = displacement vector of shaft 
^ — relative displacement 


and 


+y s ji ^ = xi +yj 


l, =( x s _x )>+(y s -y)j 

(2.91) 

n = cosOi + sin 0j 

(2.92) 

L " = L - x ) cos 9 + (y s - y) sin 6 

(2.93) 

H = C + (x s - x) cos0 + (y s - y) sin 6 = C + (x, - x) cos0 + (y 5 - y) sin 0 

(2.94) 


where: 

z = longest distance from CG to end of seal. 

2.18 Units 

The units for the English and metric systems are as follows: 

English 

• Length: in. 

• Density: lb/in. 3 

. E: lb/in. 2 

• Stiffness: lb/in., Ib/rad, in. -lb/rad 

• Damping: Ib-sec/in., Ib-sec/rad, in.-lb-sec/rad 

• Rotational speed and frequency: rad/sec 
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• Viscosity: lb-sec/in . 2 (reyns) 

• Pressure: lb/in . 2 

• Force: lb 

• Film thickness: in. 


Metric 

• Length: meters 

• Density: kg/m 3 
. E: N/m 2 

• Stiffness: N/m, N/rad, N-m/rad 

• Damping: N-s/m, N-s/rad, m-N-s/rad 

• Rotational speed and frequency: rad/s 

• Viscosity: Pa-s = N-s/m 2 

• Pressure: Pa = N/m 2 

• Force: N 

• Film thickness: microns 
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Figure 4. Program Flow Chart 
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F a = Axial force 

F f/I = Fluid film force 

F fz = Axial friction force 

5 Z = Axial seal ring displacement 

k zz = Fluid film axial stiffness coefficient 


Figure 5. Initial Equilibrium Algorithm 
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Figure 6. Spring Forces and Moments 
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Figure 11. Flow Chari of Piston Ring Wall Friction Restraining Algorithm 
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3.0 SAMPLE PROBLEMS FOR CODE DYSEAL 


The sample problems included in this section are intended to demonstrate program usage and do 
not necessarily represent seal designs. Face seal samples are included in this section. Ring seal 
sample problems are demonstrated in the verification section where output was compared against 
published data. 

3.1 Sample Problem 1: Piston Ring Face Seal Input 

The sample problem analyzes the 50-mm spiral-groove seal described in Reference 5. Key 
geometrical parameters are shown in Figure 13. The seal ring was partitioned into three elements. 
The first element was the seal ring face. The inside radius of the seal ring face was taken as the 
same as that of the other two elements, since the actual length of the inside radius of the face is 
very thin (-0.20 in.). Input is shown in Figure 14. Spring and damping coefficients were taken 
from the work reported upon in Reference 5. The shaft displacements were 0.0005 in. and 
rotations were 0.0005 radians, respectively. 

Plotted output is shown in Figures 15 through 26. In some cases, multiple plots were employed, 
such as in Figure 15, which shows the x displacement of the seal ring, XS, and the x displacement 
of the shaft, X. The shaft displacements are sinusoidal patterns, while the seal ring displacements 
are the generally lower amplitude and more irregular response, due to secondary seal friction. 

3.2 Sample Problem 2: Continuation 

Sample Problem 2 is a continuation of Sample Problem 1, for another 5 revolutions or 500 time 
steps. Continuous variables are obtained from the printed output of Sample Problem No. 1. One 
plot of minimum film thickness was made and is shown in Figure 27. 

3.3 Sample Problem 3: O-Ring Secondary Seal 

The O-ring secondary seal introduces several new variables. These include: 

• Secondary Seal Stiffness Per Unit Length which is available from O-ring catalogs. 
Figure 28 shows percent compression versus load per lineal inch of seal for various O- 
ring cross sections and durometer. This figure was extracted from a Parker O-ring 
catalog. Percent compression is obtained by calculating the squeeze of the ring 
divided by the cross sectional width of the ring. This information is also obtainable 
from the catalog. The load per lineal inch divided by the squeeze gives the stiffness 
per unit length. 
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• Secondary Seal Damping Per Unit Length. Information is not available from 
O-ring catalogs. Elastomers are light damping devices, and a reasonable number is 
1/(RSC 7c). 

• Secondary Seal Preload Per Unit Length. Generally equal to secondary seal 
stiffness per unit length, spring stiffness per unit length. 

• Radius to Secondary Seal. 

The input for this case is shown in Figure 29. 

Figure 30 shows the x displacement of the seal ring versus shaft revolutions. The displacement is 
very small - on the order of 20 x 10* 6 milli-radians. Similar results apply to the y motions shown 
on Figure 3 1 . The only exciting forces in these modes is viscous shear between the mating ring 
and seal ring, which is a small value. Figure 32 shows the axial displacement of both the runner 
and seal ring. The seal ring response is in phase and slightly magnified above the excitation. 
Rotations about the x-x and y-y axes are shown in Figures 33 and 34, respectively. The minimum 
film thickness as a function of shaft revolutions is indicated in Figure 35. Axial and rotational 
friction are shown in Figures 36 through 38. 

3.4 Ring Seal Sample Problems and Verification 

For purposes of sample problems and verification, the ring seals described by Kirk in Reference 6 were 
analyzed. Kirk used a time transient scheme that varied the fluid film forces at each time step using 
short bearing theory. That differs from the approximate analysis used in DYSEAL in which the fluid 
film is represented by constant but rotating cross coupled stiffness and damping coefficients. Kirk also 
assumed a rotor modal mass supported on springs and dashpots and determined the response of the 
modal mass. In DYSEAL, rotor motions are prescribed. 

The first case, as represented by Figure 39 (Kirk, Figure 7), included the following parameters: 

• Speed ~ 1780 rpm 

• Axial interface force = 20.9 lb 

• Ring mass = 2.13 lb 

• Clearance = 0.003 in. 

• Length = 0.904 in. 

• Viscosity = 0.8125 x 10' 7 lb-sec/in. 2 

• Friction coefficient = 0.150 

• P Ligh = 72 psig 

• Plow = 60 psig 

• Shaft radius = 2.090 in. 

• Rotor excursion = 0.0024 in. 
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The film stiffness and damping coefficients were obtained from external codes at an eccentricity ratio of 
0.5. The eccentricity was chosen on the basis of load capacity to overcome the friction forces of the 
secondary seal. The cross coupled coefficients are indicated on the input file. Figure 40. A model was 
configured that simulated the mass of the ring and provided identical wall interface and friction forces. 
The given rotor orbit was circular at an eccentricity ratio of 0.8, or a finite value of 0.0024 in. This 
corresponds to the shaft eccentricity prescribed by Kirk. The given rotor circular orbit is indicated in 
Figure 41. Orbital response of the seal ring is shown in Figure 42. There is a significant amount of 
orbit overlapping and the orbit is confined to the 0.0024-in. stimulus from the shaft. As shown at the 
bottom of Kirk's Figure 7, the orbit does overlap in a similar manner. Figure 43 shows the seal ring 
displacement as a function of shaft revolutions. Note that there is a strong half frequency component 
because of the strong cross coupling influence of the stiffness coefficients. The subsynchronous 
component further explains the orbital loops. Figure 44 shows the y displacement as a function of shaft 
revolutions. The minimum film thickness is indicated in Figure 45. From Kirk's Figure 7, the seal is 
tracking the rotor at approximately 0.5 eccentricity or with a minimum film of 1.5 mil. From Figure 
45, the median of the film variation is approximately 1.5 mil. Figures 46 and 47 show the friction 
forces in the x and y direction, respectively. The comparative results are excellent, especially 
considering the differences in problem solution methods. 

A similar analysis was conducted for Kirk's Figure 8 problem (see Figure 48). This seal was 
identical to Kirk's Figure 7 problem except that the length of the seal was reduced from 0.904 to 
0.600 in. This required the development of a new set of stiffness and damping coefficients. The 
input for this case is shown in Figure 49. The seal ring orbital response is shown in Figure 50. It 
is a complex pattern of interior looping. The minimum film thickness predicted by DYSEAL is 
shown in Figure 5 1 . Kirk indicates the seal ring tracks at an eccentricity of 0.75 or at a minimum 
film thickness of 0.75 mils, which is verified on Figure 51. The minimum film thickness on Figure 
51, however, is diminishing with revolutions and may eventually fail. Examination of Kirk's 
orbital plots reveal that the orbit is continuing to expand after the three revolutions that Kirk 
examined, and the orbit is not confined. If Kirk increased the number of revolutions, he may have 
come to the same conclusion as DYSEAL - that this ring may eventually fail by contact. 
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Figure 13. Geometry for Sample Problem 1 
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Sample Problem 1 


DYNAMIC ANALYSIS OF SEALS 


XSDISP. 
X- DISPL. 



REVS. 

95TR34-V5 

Figure 15 x Displacement versus Shaft Revolutions 
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Figure 16. y Displacement versus Shaft Revolutions 
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95TR34-V5 

Figure 17. z Displacement versus Shaft Revolutions 
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Figure 18. Film Thickness versus Shaft Revolutions 
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DYNAMIC ANALYSIS OF SEALS 


IlMIN . 
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Figure 19. Minimum Film Thickness versus Shaft Revolutions 
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Figure 20. Rotational Displacement About x Axis versus Shaft Revolutions 
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Figure 21. Rotational Displacement About y Axis versus Shaft Revolutions 



Figure 22. x Friction versus Shaft Revolutions 
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Figure 23. y Friction versus Shaft Revolutions 
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Figure 24. z Friction versus Shaft Revolutions 
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Figure 25. Friction Moment About x Axis versus Shaft Revolutions 



Figure 26. Friction Moment About y Axis versus Shaft Revolutions 
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Figure 27. Sample Problem 2 Minimum Film Thickness versus Shaft Revolutions 
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Figure 28. Typical O-Ring Data for Computing Stiffness and Preload Per Unit Length 
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Figure 29. O-Ring Sample Problem Input 



Figure 30. O-Ring Sample Problem x Displacement versus Shaft Revolutions 



Figure 3 /. O-Ring Sample Problem y Displacement versus Shaft Revolutions 
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Figure 32. 0-Ring Sample Problem Axial Displacement versus Shaft Revolutions 
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Figure 33. O-Ring Sample Problem Rotation About x Axis versus Shaft Revolutions 
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Figure 3 4. O-Ring Sample Problem Rotation About y Axis versus Shaft Revolutions 



Figure 35. ()~Ring Sample Problem Minimum Film Thickness versus Shaft Revolutions 


NAS A/CR— 2004-2 1 3 1 99/VOL5 


51 







95TR34-V5 

Figure 36. O-Ring Sample Problem Axial Friction versus Shaft Revolutions 



Figure 37. O-Ring Sample Problem Rotational Friction About x Axis versus Shaft Revolutions 
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Figure 38. O-Ring Sample Problem Rotational Friction About y Axis versus Shaft Revolutions 
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= Modal mass of rotor shaft (lb ) 

nr 

= Rotor shaft speed (rpm) 

= Horizontal or vertical indicator (1 = horizontal; 0 = vertical) 

= Rotor shaft modal mass eccentricity relative to seal clearance (DIM) 
= Seal ring mass (lb m ) 

= Seal journal radius (in.) 

= Seal radial clearance (in.) 

= Seal axial length (in.) 

= Absolute viscosity of sealing fluid (lb-sec/in. 2 ) 

= Face friction factor (DIM) 

= Seal high pressure (lb/in. 2 ) 

= Seal low pressure (lb/in. 2 ) 

= Seal liquid cavitation pressure (lb/in. 2 ) 

= Time step for transient simulation (rad) 

= Maximum time for each segment of response (rad) 

= Axial grid points for pressure profile 
= Pressure profile grid points around circumference of seal 
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Figure 39. Pump Seal Transient with Three Cycles of Motion Showing Seal 
Tracking Rotor at 0.5 Eccentricity (N = 1 780 rpm = 29.7 Hz) 
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Figure 40. Kirk’s Figure 7 DYSEAL 




Figure 41. Kirk’s Figure 7 Rotor Orbit (mil) 



Figure 42. Kirk’s Figure 7 DYSEAL Seal Ring Orbit (mil) 
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Figure 43. Kirk’s Figure 7 DYSEAL x Displacement (mil) 



Figure 44. Kirk ’s Figure 7 DYSEAL y Displacement (mil) 
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Figure 45. Kirk’s Figure 7 DYSEAL Minimum Film Thickness (mil) 


DYNAMIC ANALYSIS OF SEALS 
CHECK AGAINST G. KIRKS RESULTS, FIG. 7 
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Figure 46. Kirk’s Figure 7 DYSEAL x Friction (lb) 
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Figure 47. Kirk ’s Figure 7 DYSEAL y Friction (lb) 
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Figure 48. Pump Seal Transient for a Reduced-Length Seal Showing Seal Ring 
Tracking Rotor at an Eccentricity of z - 0.75 (N - 1 780 rpm = 29. 7 Hz) 
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Figure 49. Kirk’s Figure 8 DYSEAL Input 







4.0 VERIFICATION FOR CODE DYSEAL 


Several methods of verification were accomplished. Internal checks were made against closed- 
form solutions. Some mass, spring, and damper vibration problems were examined and compared 
against closed-form solutions. Also, comparisons were made against experimental data available 
in the literature. Ring seal verification was presented in the preceding section. 

4.1 Internal Checks 

An example of internal checks of the code was evaluation of spring forces and moments. The 
numerical approach used in the code can be compared against a closed-form solution. For the ith 
spring, the moment is 

=r S pXr‘ =(z !p -z cg )k + (R sp cos0 i t + R sp sin0,j)x(-k s z i k) (4.1) 


where: 


= moment from ith spring located at 
F s ! = force from ith spring 
z sp = z distance to ith spring 
R^ p = spring radius 

= angular coordinate to ith spring 
Zj = displacement of ith spring in z direction 
z, 0 = z distance to CG 

'-g 

The seal ring motion in the z direction, z i3 is given by 

z, =z s +R sp (p s sin9 1 -a s cose i ) (4.2) 

where the variables have been previously defined. 

After substituting Equation (4.2) into Equation (4.1), expanding and summing over all springs, the 
following equation results: 

M s = (R sp k s z s loose, + R‘ p k,p s I sin 9 ; cosO, - R s 2 p k s a s E cos 2 9, )] 

-(R;k s z s Esin 6j + R 2 p kJ3, Esin 2 6, -R 2 p k s a s Esin9, cos9,)i 

where the j component equals the moment about the y axis and the i component equals the 
moment about the x axis. 

Consider six springs with the first spring starting at an angle of 10°. The summation coefficients 
are indicated on Table 3. 
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Table 3. Summation Coefficients 


Kjp = 2; N sp = 6; k = 100 lb/in. 


0 

sin0 

COS0 

sin0 cos0 

sin 2 

cos 2 

10 

0.1736 

0.9848 


0.0301 

0.9698 

70 

0.9397 

0.342 


0.883 

0.117 

130 

0.7660 

-0.6428 


0.5866 

0.4132 

190 

-0.1736 

-0.9848 


-0.301 

0.9698 

250 

-0.9397 

-0.342 


0.883 

0.117 

310 

-0.7660 

0.64281 


0.5868 

0.4132 

2 

0 

0 

0 

3.0 

3.0 
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Referring to Equation (4.3), 

K x „ = R, p k s Xcos0 = 0 

K pc =R" p k s Isin0cos0 = O 

K aa = R«pk s Icos 2 0 = 2 2 (100)(3) = 1200 lb / in. 

Similarly, 

K z p = K a p = 0, Kpp - 1200 

Also, the axial stiffness - 6x 100 = 600 lb/in. The spring stiffness matrix as computed by the 
program is as follows: 

paste up little table here 


This checks precisely with the closed-form solutions. Similar results were obtained for varying 
the number of springs and the independent stiffness values. 

4.2 Mass, Spring, Damper Vibrations 

Several mass, spring, and damper vibration problems can be used to check out portions of the 
code. First, consider a forced vibration problem, as depicted in Figure 52. The base represents 
the shaft; the mass represents the seal ring. The base is the source of excitation, and the response 
of the seal ring or mass, M, is desired. The initial parameters tested were: 

k = 10,000 Ib/in. 
t — 10 lb-sec/in. 

M = 10 lb = 10/386.4 = 0.02588 lb-sec 2 /in. 

X = 0.002 in. 
co = 1000 rad/sec 
C = 10 lb-sec/in. 

As derived from Thomson [7], the maximum relative displacement, Z, (= y - x) is given by 

mco 2 X (0.02588)(1000) 2 (0.002) 

-J(k-mar) 2 +(Cco) 2 ^(l 0, 000 - 0. 025 9( 1 000) 2 f + ( 1 0(l 000)) 2 (4.4) 

= 0. 002758 in. = 2. 758 mil 
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Z represents the maximum difference between the amplitude of the mass, y, and the excitation, X. 
As shown on the computer output graph (Figure 53), the measured difference equals 2.74 mil. 
The phase angle as a function of frequency and damping is shown on Figure 54. 


03 _ co _ 1000 

r.A (~k. / 10,000 

n "V m V 0.0258 


C C 
C c 2mco„ 


10 

2(0.02588)(622) 


= 0.3106 


(4.5) 

(4.5) 


From Figure 54, the phase angle, c|), is approximately 135°. The computed value is estimated to 
be 136° as measured from the output curves of Figure 53. Considering graphical interpretations, 
the corroboration is excellent. Similar results were obtained using the ring seal option of the code 
exciting the shaft in the x direction, as shown on Figure 55. 

4.3 Verification Against Data in the Literature 

Di Russso [8] did extensive dynamic testing of spiral-groove face seals. The seals were subjected 
to constant load rotation, with installation runout of the seal seat (rotating member) and then to 
constant load rotation with installation runout plus an axial excitation of 50 fim (2 mil) at 100-Hz 
frequency. Figure 56 schematically shows the seal seat vibrational modes. The installation 
misalignment of the seal seat was approximately 35 p,rad about both the x and y axes. Tests were 
run with and without secondary seals. 

Figure 57 shows response of the seal at 14,000 rpm without a secondary seal in place. The film 
thickness frequency is approximately 6 times synchronous for the 14,000-rpm case. The case was 
simulated by determining the stiffness and damping characteristics of the spiral groove and 
establishing the physical characteristics of the seal ring. The input for the DYSEAL run is shown 
in Figure 58. The minimum film thickness in mils versus shaft revolutions, computed by the code 
is shown in Figure 59. The 6 per rev frequency is shown in Figure 60. The results are nearly 
identical to the steady mode without the axial excitation. The implication is that the seal ring 
tracks the exciting shaft perfectly. The input for the computer studies is identical to Figure 58, 
except that Z0 is give a value of 0.002. Computer results of film thickness are shown in Figure 
61. The film thickness shows a definite trace of the excitation frequency. A blown-up view of the 
film thickness is shown in Figure 62, and the six times synchronous frequency is clearly 
discernible. The axial displacement of the seal seat (ZS) and seal ring (Z) are shown in Figure 63. 
They are in unison, confirming the tracking capability of the seal as experienced on test. The 
variations indicated by the film thickness curves do not show in the axial mode but are indicated 
by the rotational response about the x and y axes. Figure 64 shows the rotational response about 
the x axis superimposed on the pure sinusoidal excitation. The jagged sine curve provides the 
differences between excitation and response. In general, the computational results agree very 
favorably with the experimental data. 
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Figure 52. Mass , Spring, and Damper System 
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Figure 53. Single-Dcgree-of- Freedom Forced Vibration 
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Figure 54. Phase Angle as a Function of Damping and Frequency 



95TR34-V5 


Figure 55 . Ring Seal Option: Single-Degree-of-Freedom Forced Vibration Problem 
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a) Steady Seal Seat Mode with True Seal Seat 



b) Sinusoidal Seal Seat Mode with True Seal Seat 
Figure 56. Schematic Showing Seal Seat Vibrational Modes 
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Figure 57. Film Thickness as a Function of Time (Probe 1) for Inward-Pumping 
Spiral-Groove Seal (No Secondary Seal) and Steady Seal Seat Mode 
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Figure 58. Input for Spiral-Groove Seal; 14,000 rpm, No Axial Excitation 
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Figure 59. Results of DYSEAL Analysis; Film Thickness versus 
Revolutions (Steady Seal Seat Mode) 
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a) Steady Seal Seat Mode 



b) Sinusoidal Seal Seat Mode; Amplitude = 50 pm (2 mil); 
Frequency = 100 Hz 
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Figure 60. Film Thickness; Sinusoidal Axial Vibration 
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Figure 61. DYSEAL Film Thickness; Sinusoidal Axial Vibration 
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Figure 62. DYSEAL Magnified View of Film Thickness; Sinusoidal Axial Vibration 
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Figure 63. Axial Motion of Shaft and Seal 



Figure 64. Rotational Response About x Axis for Axial Sinusoidal Excitation 
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5.0 DESIGN MODEL DESCRIPTION OF CODE KTK 


The labyrinth seal design model is an expansion of the knife-to-knife (KTK) analyses reported in 
the literature. In such approaches, one-dimensional flow parameters in the knife throats are 
computed and linked together by a total pressure loss calculation. Flow coefficients are used for 
individual knives or groups of knives to account for the vena contracta in knife throats. Velocity 
head carry-over from upstream knives is accounted for by reducing the head loss between knives 
based on the flow expansion angle. 

The design model is similar to previous KTK approaches except that the loss for each knife is 
broken down further into three losses as shown in Figure 65: 

• Contraction: stations 1 to 2 and 4 to 5 

• Venturi and friction: stations 2 to 3 and 5 to 6 

• Partial or full expansion: stations 3 to 4 and 6 to 7. 

The three loss coefficients aid in understanding the effects that each seal parameter has on the 
pressure drop across a knife. They reflect the types of pressure drops that the flow experiences 
more specifically than a knife flow coefficient used in previous KTK models. Specifically, the 
design model consists of : 

• One-dimensional flow calculation at three locations for each knife 

• Calculation of the individual loss coefficient values from the flow and geometric 
conditions 

• Correction of loss coefficient values due to the presence of adjacent knives 

• Linking of the total pressure between stations based on the total pressure drop from 
the loss coefficient and local velocity head. 

Figure 66 lists the basic flow equations used in the design model. 

Early in the development of the design model, knife loss coefficients were specified in the input 
data. The loss coefficients were then corrected by applying the model to the data base. The 
corrected coefficients were correlated against the independent seal parameters using a linear 
regression analysis. This iterative process was continued until equations were obtained which, not 
only, fit the data but were also physically relevant. The correlating equations were put into the 
model code and then the model accuracy was assessed by comparing the calculated flow results 
with the flow results in the entire data base. 

Detailed information about the design model is summarized in the following paragraphs. This 
includes a discussion of the parameters considered and the model development for single knives, 
straight seals and stepped seals. 
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5.1 Parameters Considered 


Based on the literature survey and previous experience at Allison [9], the parameters selected for 
consideration in the design model were those listed in Table 4. Schematics are given in Figures 67 
and 68 which define these parameters. The effect of rotational speed was not included because of 
the lack of accurate and complete experimental data for labyrinth seals at knife tip speeds typical 
of those found in gas turbine engines. 

The parameters listed in Table 4 were incorporated into the model to define the seal geometry. 

The geometric parameters were combined to obtain the data correlations in a manner that would 
make the correlations physically relevant. The ranges of these geometric parameter combinations 
in the data base are given in Table 5. The design model is valid primarily within these ranges. 

5.2 Single Knife 

Analyzing the flow data for a single knife affords the advantage of isolating the individual knife 
losses from the influence of adjacent knives. For the design model approach the three loss 
coefficients for the knife were separated using overall flow characteristic data. This was done by 
first noting that the expansion loss (IQ would be unity because of the complete expansion 
downstream of a single knife. The other two loss coefficients, contraction (Kc) and venturi with 
wall friction (K v f), were separated by considering data from several authors who varied the 
leading-edge sharpness and the thickness of the knife tip independently. It was assumed that the 
leading-edge radius affected only Kc and that the knife tip thickness affected only Kvf. The 
resulting correlations for Kc and Kvf are given in Figures 69 and 70. In these correlations, the 
local flow factor, <J>, is an independent parameter. Several values for knife tip radius (KR) are 
listed in Figure 69 from the five data sources. 

Single knife flow data for slanted knives were also analyzed to determine the effect of knife angle. 
Figure 69 gives the results. For K9 < 90°, it was found that a parametric function of K0 added to 
Kc for a 90° knife gave the best fit of the data. For K0 > 90°, i.e., K(3 = 0, data from Idel’chik 
[14] indicated that a multiplier applied to Kc for a 90° knife would be the most appropriate. 

The relationships for the three loss coefficients as summarized in Figure 69 represent one 
approach that could be taken to correlate overall flow characteristics of single knife data. In 
combination, they accurately define the overall flow characteristics for the data sources and 
parameter ranges considered. The individual correlations involve the interactions of seal 
geometry with the flow parameters in a physically realistic way so that the overall results are valid 
for interpolated and some limited extrapolated data. 

5.3 Straight Seals 

Multiknife seals are analyzed in the design model by linking the triplet losses for each knife in 
series so that a seal has a total of 3 times KN losses. The overall pressure loss is a summation of 
the individual total pressure losses. The losses are calculated sequentially starting with the known 
inlet pressure because the loss coefficients and Mach number are functions of the local total 
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pressure through the parameter <\>. For a straight seal, there is a carry-over of the velocity head 
from one knife to the next. This carry-over affects the Kyf and fQ of the upstream knife and the 
Kc and K v f of the downstream knife. Thus, all the loss coefficients of a multiknife straight seal are 
influenced except the Kc of the first knife and the K« of the last knife. The approach followed in 
multiknife seals is to determine the three loss coefficients for a given knife location from single 
knife correlations (Figure 69) and then correct them for the influence of adjacent knives. The 
correction is based on the characteristics of the expansion angle of the jet discharging from the 
clearance gap over a knife. This approach has been discussed by Abramovich [16] and utilized by 
Komotori and Miyake [l 1] in their KTK model. Figure 71 shows a schematic of a straight seal 
with the expansion angle, a, defined. The flow expands until it impinges on the front edge of the 
next knife. The maximum downstream flow height is CL + 6 so that the expansion area ratio is 
(CL + 5)/CL instead of (CL + KH)/CL if the next knife were not present. This jet expansion ratio 
not only represents the amount the flow expands from the upstream knife but also the contraction 
into the downstream knife gap. The equations for 6 in terms of a and the other geometric 
parameters are given in Figure 71. 

To incorporate the effect of a on the three loss coefficients, relationships given in Figure 72 were 
used which are recommended in the literature [17] for various types of losses. The ratios Ao/Ai 
and A 0 /A 2 are simply the ratio CL/(CL + 6) relative to the upstream and downstream sides of a 
given knife, respectively. 

The expansion angle, in general, will vary from knife to knife because the pressure ratio varies. 
This was observed in the flow visualization test results. This expansion angle variation, however, 
was not modelled because of the lack of good, complete seal data with interknife pressure data. 
Future design model development could be done to include variation through the seal based on 
results from analysis model calculations and/or test data. 

The equations in Figures 71 and 72 were incorporated into the design model with a as an input 
value. Results were then obtained for the straight seal geometries in the data base using a range 
of a values. Comparing model results with the test performance data yielded the correct a 
values. Table 6 summarizes the range of values obtained from the various data sources. The a 
range obtained from Komotori’ s data [11] compares well with the value of 6° reported in a 
discussion of their paper. A linear regression analysis was performed on the a results. The 
equation obtained in given in Figure 73. 

The effect of land roughness was included in the model by adding a frictional head loss term to 
the venturi with a wall friction loss coefficient (Kf) in the model for a smooth land (see Figure 73). 
Also, the flow area in each knife throat was increased to account for the increase in clearance due 
to the land roughness. The frictional head loss coefficient was determined from published rough- 
pipe turbulent flow friction factor correlations and a length equal to the knife pitch for knives 2, 3, 
and equal to knife-tip thickness for the first knives. 


NAS A/CR— 2004-2 1 3 1 99/VOL5 


77 



5.4 Stepped Seals 


Stepped seals are designed to minimize carry-over from one knife to the next. Accordingly, one 
model approach could be to simply derive a new correlation for the expansion angle in terms of 
cavity dimensions. This method, however, is not acceptable because stepped seals flow both 
more and less than equivalent straight seals without carry-over depending primarily on the 
clearance. An expansion angle approach can only account for a flow equal to or greater than that 
without carry-over. 

Physically, the flow between knives in a stepped seal does carry over some of the velocity head to 
the next knife. But while the intervening flow path dissipates some of the velocity head it also 
affects how the flow enters the next knife and thereby influences the loss coefficients of that knife. 
The complex flow patterns involved would make correction correlations to the loss coefficients 
difficult to determine accurately. Consequently, a simpler approach was taken to introduce an 
area correction factor (XMUL) for a knife throat downstream of a step. This factor is a multiplier 
on the flow area and can be less than or greater than unity. It accounts for carry-over, additional 
pressure loss in flowing between the knife face and step which is important for small distances to 
contact (DTC), and flow distortion into the next knife throat. 

A correlation for XMUL was obtained through a procedure similar to that followed for a for 
straight seals. Model results were calculated for stepped seal configurations in the data base for a 
range of input XMUL values. Comparing these results with test data yielded the correct XMUL 
value. The area multiplier was found to vary from 0.55 to 1.32. A correlating equation for 
XMUL in terms of the various geometric parameters was derived using a linear regression 
analysis. This was done first for STLD flow direction, backward facing stator steps, because 
there were 62 configurations for STLD flow direction compared to 15 for LTSD flow. A 
correlation for the LTSD flow direction was obtained by comparing the STLD equation to the 
LTSD data and deriving a correction expression. This approach gives the best chance to 
extrapolate the narrower parameter ranges for the LTSD data. Figure 74 gives the two 
correlating equations for XMUL and respective parameter ranges. 

Roughened surface land effects for stepped seals were handled in the model using a procedure 
similar to that used for straight seals, i.e., adding a friction head loss term to Kyf for a smooth wall 
plus increasing the throat area due to roughness. However, the length of the roughened passage 
was taken to be the knife-tip thickness to give the best agreement between the model and test 
data. 

5.5 Design Model Optimization 

The design model is an abbreviated analysis tool. Design is typically accomplished using this 
model by (1) determining the overall design constraints, (2) selecting the allowable range of each 
parameter to meet design and model constraints, (3) using the design model to calculate the 
leakage flow rate for a matrix of possible seal configurations, and (4) optimizing the seal design 
from the performance matrix, i.e., finding the seal geometry with the lowest leakage. This process 
has been automated by coupling the design model with a numerical optimization algorithm. As a 
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result, a minimum amount of input information is needed to optimize a seal configuration using 
the design model. In the following subsections, a brief description will be given for the 
parameters considered in the optimization algorithm. 

5.5.1 Optimization Parameters 

The parameters considered in the optimization code are listed in Table 7. The parameters are of 
three types: (1) input parameters which define the seal configuration but are not optimized, (2) 
optimized parameters, and (3) constrained correlation parameters. The input parameters are not 
optimized because of their nature (T, Pi, Pr) or the design defines their value (CL, KR, K9, Lmax, 
Hm ax, DTC, DIA, KPmin). The parameters and H,^ are optional and constrain the 
calculations only if input. The optimized parameters listed in Table 7 are of two types: continuous 
and discrete. These types are handled differently by the optimization algorithm. The third type of 
parameter constrains the selection of the best design so that the various correlations in the design 
model are not extrapolated. 

5.5.2 Optimization Algorithm 

Determining the best seal design is an iterative selection procedure which is the subject of a 
branch of mathematics known as optimization theory. The characteristics of the problem solved 
determine the type of theory; the most general is nonlinear constrained optimization. In this case, 
a nonlinear objective function is optimized with respect to the design requirements, also called 
independent variables or parameters, that are subject to equality or inequality constraints which 
are also nonlinear functions of the independent variables. The method selected for the design 
model optimization code involves the use of a penalty function to convert the constrained 
optimization problem into an unconstrained one which is solved using the Fletcher-Power- 
Davidon variable metric method. This algorithm includes a parabolic cubic spline fit search 
routine to locate the optimum. This approach is reliable even for erratic functions often 
encountered in design problems. 

The optimization algorithm described applies only to continuous variables, i.e., the first six 
optimized parameters listed in Table 7. The discrete variables were optimized by trial and 
comparison in which the entire matrix of these variables is considered. The code performs the 
continuous variable optimization for each set of discrete variable values and the overall optimum 
design is selected from the individual optimum designs. 

Constraints have been included in the algorithms to ensure that the optimized seal configuration 
satisfies the design requirements. The constraints imposed are given in Table 8. Constraints on 
the discrete variables (KN, seal type, and flow direction) simply limit the matrix of values 
considered in the trial and comparison procedure. Constraints on the other variables are imposed 
by adding inequality penalty functions to the functions being optimized. A penalty function equals 
zero if the design meets a given constraint . It is greater than zero if the constraint is violated and 
the penalty varies parabolically with the magnitude of the violation. Each continuous variable 
constraint has one penalty function associated with it. 
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Compressible flow equations applied at the knife throat are: 
the isentropic pressure relationship 


Y 



combined with the compressible flow equation of SAINT VENANT-WANTZEl 



in the form 
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Figure 66. Basic Flow Equations Used in the Design Model 
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Figure 6 7. Seal Nomenclature for Straight Seals 
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Figure 68. Seal Nomenclature for Stepped Seals 
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Table 4. Parameters in the Design Model 

Geometric parameters for straight and stepped seals 

o Knife height (KH) 
o Knife pitch (KP) 
o Number of knives (KN) 
o Knife angle (KG) 
o Knife tip thickness (KT) 
o Knife taper angle (KB) 
o Knife tip leading edge radius (KR) 
o Clearance (CL) 
o Surface roughness (e) 

Additional parameters considered for stepped seals 

o Step height (SH) 
o Distance to contact (DTC) 
o Flow direction (USD or STLD) 

Flow parameters 

o Overall pressure ratio (Pr) 
o Inlet stagnation pressure (Py) 
o Fluid temperature distribution (T) 
o Flow rate (w ) 

95TR34-V5 
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Table 5. Parameter Ranges of Data in Labyrinth Seal Data Base 


Pa ramete r 


Sea 1 

type 

Stepped 

sea 1 

Single 

knife 

Straight 

seal 

STLD dir. 

LTSD dir 

KN 

mi n 

1 

2 

2 

2 


max 

1 

12 

6 

6 

KT/CL 

min 


0.21 

0.21 

0.50 


max 

3.3 

4.4 

2.64 

1 .50 

KG 

mi n 

30 

60 

50 

50 


max 

90 

90 

90 

90 

KH/Cl 

mi n 

- 

2.7 

5.1 

5.1 


max 

- 

31 .3 

29.4 

28.0 

KP/CL 

mi n 

- 

4.0 

6.4 

9.2 


max 

- 

56.3 

53 

40 

c/( 2CL) 

mi n 

0 

0 

0 

0 


max 

0 

0.030 

0 

0.030 

SH/CL 

mi n 

- 

- 

2.0 

4.0 


max 

- 

- 

29.4 

12.5 

DTC/CL 

mi n 

- 

- 

0.85 

4.1 


max 

- 

- 

40 

19.4 

( KP-KT ) /CL 

mi n 

- 

3.5 

6.2 

8.9 


ma x 

- 

55.0 

51 .8 

38 . 5 
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K -1.0 

e 

K t - f (4>, KT/CL) 
vf 


where 


Figure 70: 0.77 < KT/CL < 3.3, but good for 0.0 < KT/CL < 3.3 

[Derived from Kearton and Keh data for K^ = 0.70 (KR very 


smai i ) ] 



KR - i n 
0.0 


0.00167 

0.005 

0.005 

0.010 


EXP 



C 2 ♦ 


2/ CL 


0.25 


from 



KEARTON & KEH [10] 


Allison 

K0M0T0RI & MIYAKE [I I] 
HARRISON [12] 

CAUNCE & EVERITT [13] 


K 

c 


K 

c 


K c @ 90° 

K c <a 90" X [1 ‘ " C 4 <K6 “ 90 ° )] 


S @ 90° + S O. - SIN (K9) ] 

[from Allison plus MEYER AND LOWRIE 


for K9 = 90° 
for K0 >90° 

[from IDEL'CHIK [14] ] 
for 30° < K0 < 90° 
data [15] ] 


NOTE: Ke is actual front surface angle relative to the flow direction so that 
Ke = 90° + KB/ 2 when the specified knife angle is vertical or beyond, 
KO > 90°. 

C = constant 

n 
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Figure 69. Loss Coefficient Correlations for Single-Knife Seal 
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4 - ( KP-KT) /[Tan kB + (1/Tan a)] KG = 90 deg 

(KP-KT) Tan a a < KG < 90 deg 

95TR34-V5 

Figure 71. Schematic of the Flow Expansion Angle for a Straight Seal 
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Figure 72. Effect of Upstream and Downstream Area on Loss Coefficient 
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Table 6. Expansion Angle (a) Determined by Correlation 


6-8 deg 
4-6 deg 
2-4 deg 
4-5 deg 
2-4 deg 


JET EXPANSION ANGt F 


o Caunce and Everett, 6 knife 
o Komotori 2, 4, 8, and 10 knife 
o AGT 4 knife 
o AGT 8 knife 
o AGT 4, 5 knife slanted 



for 


0.54 


< ( KP-KT) 
KH 


< 4.0 


Average deviation = 25% 

C 6 = constant. The value of this constant is given in the program 
listing for the Design Model, Appendix M D". 


WALL ROUGHNESS 


K v f " K vf smooth (Correction for upstream and downstream knives) + K 


f rough 


where 


K f rough - f < C/H ’ Re - KP > 


A = a ! ^ + - 

t t smooth V CL 
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Figure 73. Straight Seal Correlations in the Design Model 
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STEPPED SEAL AREA MULTIPLIER. XHUL 


STEP Flow Direction 


c c c 

XMUl = C 7 ( OTC/CL) (KT/CL) ( PTC / ( KP-KT) ) 9 (KH/CL) 10 
•" ( (KP-KT) /KH) 11 (SH/CL) Cl2 /^(dTC/CL) 2 7”c _ 

0.85 < OTC/CL < 40, 0.21 < KT/CL < 2.6, 0.09 < DTC/(KP-KT) < 1.0, 

5.1 < KH/CL < 19.4, 1.16 < (KP-KT) /KH < 1.76, 2.0 < SH/CL < 29.4 

LTSD Flow Direction 

XMUL = XMUL stLD C 14 ^ KH/CL ) 

4.0 < DTC/CL < 19.4, 0.50 < KT/CL < 1.5, 0.35 < DTC/(KP-KT) < 0.50 

5.1 < KH/CL < 28, 1.02 < (KP-KT)ZKH < 1.9, 4.0 < SH/CL < 12.5 

Note. The limits on the seal parameters result from the range of 
geometries used in developing the correlation equations. 


WALL ROUGHNESS 


K r ~ K ' K 

vi vf f rough 


K 


rough 


f(c/H, Re, KT) 


A = A 


t 


smooth 



Figure 74. Stepped Seal Correlations in the Design Model 


the seal 
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7 able 7. Design Model Optimization Parameters 


Input Parameters 

Straight and Stepped Seals 

Clearance (CL) 

Temperature (T) 

Inlet total pressure (P ) 

Pressure ratio (P ) 

R 

Knife radius (KR) 

Knife taper angle (KB) 

Maximum axial lenqth {L )* 

max 

Stepped Seals Only 

Maximum seal height (H )* 
Distance to contact (DTC) 
Maximum or minimum diameter (D 
Minimum knife pitch (KP 


Optimized Parameters 

Continuous Variables 

Knife height (KH) 

Knife pitch (KP) 

Knife tip thickness (KT) 

Knife angle (Ke) 

Roughness (e) 

Step height (SH)** 

Discrete Variables 

Seal type (straight, stepped) 
Number of knives (KN) 

Flow direction (USD, STLD ) ** 


, D . ) 
max min 

mi n^ 

(= 2X maximum allowable axial travel) 


* Optional 
** Stepped seals only 


Constraining Correlation Parameters 


Straight Seals 


Stepped Seals 


KT/CL 

Ke 

( KP-KT ) /KH 
(c - 30) /CL 


KT/CL 

Ke 

( KP-KT) /KH 
DTC/CL 
SH/CL 
KH/CL 

(c - 30) /CL 
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Table 8. Constraints Imposed in Design Model Optimization Code 


Overall length (L )* 
max 

Overall height (H )* 
max 

Minimum and maximum limits for all optimized parameters, i.e., 

KH, KP, KT , K6, c, SH, KN**, seal type**, flow direction** 

Minimum and maximum limits for correlation parameters to avoid 
extrapolation beyond data ranges in the data base, i.e., 

KT/CL, K0, (KP-KTJ/KH, DTC/CL, SH/CL, (c - 30)/CL 
* Optional constraints 

** Discrete variable constraints not imposed by inequality penalty functions 

95TR34-V5 
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6.0 COMPUTER PROGRAM FEATURES OF CODE KTK 


The design model computer program described in this report has several features to enhance its 
use. The program can be run in an analysis mode by defining the labyrinth seal geometry to be 
considered and executing the program. Also, the program can be run in an optimization mode by 
defining certain parameters and letting the program determine the optimum values for the 
remaining ones. The first mode will be referred to as the design model code and the second one 
as the design model optimization code. Various features of these two parts will be described 
separately in the following paragraphs. 

6.1 Design Model Code Features 

Features available in the design model code include: 

• The input is abbreviated where possible. 

• An override is available for many of the loss coefficient parameters. 

• A function loss is available instead of, or in addition to, the three loss coefficients. 

• Straight seals, stepped seals, or a mixed combination of the two can be analyzed. 
(Steps can be for either increasing or decreasing diameter.) 

• Geometric parameter values can be varied from knife to knife. 

• Two-dimensional or three-dimensional seals can be analyzed. (Calculations for two- 
dimensional seals are important to compare model results with nonrotating labyrinth 
seal rigs that utilize a rectangular test section). 

• Various calculation options can be selected. The program can calculate: 

- Seal pressure distribution for a given flow rate 

— Seal pressure distribution and flow rate for a given overall pressure ratio 

- A flow characteristics curve (<{> versus Pr) 

• A flow characteristic curve can be punched out for input to plotting routines or to 
flow network solution codes. 

• A summary is printed out of the various parameters with their applicable ranges. (If a 
parameter is outside its range, a warning message is printed, and the calculations are 
continued with either the input parameter value or one at the end of the range 
depending on whether or not extrapolation is considered acceptable for that particular 
parameter.) 

6.2 Optimization Code Features 

The optimization code is the design model code coupled with a driver routine. The latter 
calculates the independent parameter values to be used in the design model to search for an 
optimum configuration. Features of the code are: 

• Constant geometry straight and stepped seals can be considered. However, variable 
parameters from knife-to-knife or mixed straight and stepped seal geometries cannot 
be considered. 


NAS A/CR— 2004-2 1 3 1 99/VOL5 


93 



• Each independent parameter has a default range which may be overridden. Even the 
correlation parameter ranges may be overridden if desired. 

• An independent parameter may be held constant (by inputting both its minimum and 
maximum values equal to the one desired). 

• An optimum configuration may be determined for both seal types and both directions 
for the stepped seals. Any subset of these may be considered. 

• Before optimization is attempted, the parameter values and ranges are checked to be 
sure a solution is possible, e.g., a solution is impossible if Lmax is less than the minimum 
KP divided by the maximum KN. If a solution does not exist, information is printed 
describing the problem and the execution of the data set is halted. 

• Intermediate output information is given for each combination of discrete variables 
employed. This output information includes algorithm parameter values, derivatives of 
the optimized function with respect to each continuous variable, and comparisons of 
the continuous variable values with the allowable ranges. 

• Final output information includes sensitivity results for each discrete variable step and 
summary data of the optimum seal configuration designated/ 

The output information not only defines the optimum seal configuration but indicates the 
effect, if any, of imposing each constraint. Also, the improvement in decreased leakage of 
the optimum configuration compared to the other possible configuration is given. This 
information can be used to assess the penalty caused by each limiting constraint and the 
penalty for choosing an alternate design. 

6.3 Description of Output 

The design model code printed output provides a detailed description of the seal geometric 
parameters and the predicted aerodynamic performance. A listing of a sample output data set is 
given in Appendix A. The output associated with a seal that is being optimized differs 
significantly from a nonoptimized seal. Therefore, the two outputs will be discussed separately. 

6.3.1 Nonoptimized Output 

A description of the output corresponding to sample data sets 1 through 3 in Appendix A is 
presented in the following paragraphs. An overview is presented first, followed by a more 
detailed description of the output. 

The first section of the output echoes some of the parameters input. The second section, “KNIFE 
GEOMETRY DATA,” lists the geometric parameters associated with each knife of the seal. The 
third section, labelled “INPUT DATA RANGE CHECK,” records the results from a check of the 
input data against the data ranges used in the design model correlations. Warnings are issued 
when input data forces an extrapolation outside the empirical data range. If a seal parameter is 
outside the empirical range, this output section is printed before the first section described above. 
The next section of the output lists the aerodynamic parameters for each of the three flow 
“stations’ associated with each knife. The values shown correspond to the choked flow condition 
unless a seal flow or pressure ratio is specified in the input. The fifth and final output section is 
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labeled “FLOW CURVE” and prints the values which make up the flow curve for the seal (if 
applicable). The flow curve (a function of pressure ratio and WRT/PTA) is printed out in 
standard and elliptical coordinates. 

The “KNIFE GEOMETRY DATA” is printed for the three stations of each knife. The first 16 column 
headings correlate closely with the input variable names. The remaining column headings, beginning 
with “DEL C,” identify the values for other parameters determined by the code as listed below: 

• DEL C : Radial expansion of flow jet between knives used to calculate contraction area 
ratio. 

• DEL E: Radial expansion of flow jet between knives used to calculate expansion area 
ratio. 

• AREA MULT: Area multiplier on flow area for stepped seals. 

• ALPHA: Jet expansion angle used for carry-over calculation. 

The “INPUT DATA RANGE CHECK” is printed after the “KNIFE GEOMETRY DATA” 
except when a knife parameter is outside the empirical data range. In such cases this section is 
printed at the beginning of the output to draw attention to the warning messages printed. The 
output consists of the knife parameter being considered, the minimum value of the empirical data, 
the value input (or calculated) for the knife, and the maximum value of the empirical data. 
Asterisks are printed near the minimum or maximum value of the knife parameter when the 
minimum or maximum value (respectively) of the data base has been exceeded. A warning 
message is printed to this effect. 

The fourth section of output lists the aerodynamic performance data predicted by the design model for 
each knife station. The column headings are defined below beginning with the third column. 

• AREA: The flow area over the knife tip calculated by the equation: 

A = n * DIA * CL * AREA MULT (in. 2 ). 

• TEMP: The temperature of the air entering the seal as input (°R) 

• WRT/PTA: Flow function W * VTin / Ptin * Aref(lbmVR / sec/ lbf) 

• WRT/PSA: Flow function W * VTi n / Psin* Aref(lbmVR / sec/ lbf) 

• 4FL/D: Friction factor calculated for the seal land 

• KFACT: Corrected dynamic head loss factor (K-factor) 

• LOSS TYPE: Pressure loss mechanism (sudden contraction, “long hole,” sudden 
expansion) 

• PT: Total pressure (lbfiin.**2) 

• PS: Static pressure (lbf/in.**2) 

• MN: mach number. 

• KFACT METHOD: Method of calculating total pressure loss (total pressure minus 
static pressure of 1/2 * density * velocity* *2) 

• PARM: Correction parameter 

• AMUL: Multiplier on the “XKUNC” to calculate corrected dynamic head loss factor 
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• ADDER; Quantity added to the quantity “XKUNC” times “AMUL” to calculate 
corrected dynamic head loss factor 

• XKUNC: Uncorrected dynamic head loss factor 

• MOD AREA: Modified flow area calculated by multiplying “AREA MULT” (from 
KNIFE GEOMETRY DATA printed output) by “AREA.” 

The final section of printed output is labelled “FLOW CURVE” and prints the coordinates of a 
curve relating W^/Tin / PtAref to seal pressure ratio. Also shown are coordinates for curves 
relating a function of pressure ratio to functions of W^Tin / PtAref in elliptical coordinates. 
Coordinates for a flow curve relating WjTin / PtAref and pressure ratio are also “punched” out 
on Fortran I/O unit 7. The row headings are defined below: 

• PR: Seal pressure ratio 

• PHI; Flow function (Wy/Tin / PtAref ) 

• 1 1 ./PR* *2: One minus the reciprocal of seal pressure ratio squared 

• PHI* *2: Flow function (W^Tin / PtAref ) squared 

• R/G * PHI* *2: Gas constant for air divided by the gravitational constant times flow 
function W^ffin / PtAref ) squared. 

6.3.2 Optimized Output 

A description of the output corresponding to sample data set 4 in Appendix A is presented in the 
following paragraph. 

The first section of output prints the input values for the variables on input record type 7. The 
second section of the printout lists some of the information input on record types 1 through 4. 

The third section, “KNIFE GEOMETRY DATA,” lists the geometric parameters associated with 
each knife of the seal for the maximum number of knives allowed by input. The fourth section 
lists some information relating to the first optimization iteration. The fifth output section 
(“KNIFE GEOMETRY DATA”) lists knife geometric data for the initial optimization 
configuration. The resulting calculated flow rate is then printed followed by a table of 
aerodynamic parameters for each knife station for the configuration being considered. The 
seventh section lists the optimum value for each seal parameter with the range associated with 
each (set by input or default). Parameters that are “binding constraints” are flagged. Output 
sections 4 through 7 are repeated each time the number of seal knives is indexed. The last section 
of output lists the seal input parameters that are not optimized and summarizes the optimized 
knife configurations. The optimum seal configuration (i.e., minimum leakage flow) is flagged. 
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APPENDIXA 
CODE OUTPUT 
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Mixed Labyrinth Seal 
Data set #2 
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